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ABSTRACT 



Cin Context. In a previous paper (Paper I) we showed that multi-phase gaseous halos of late-type spiral galaxies, detected in the radio continuum, 

<»>>, in Hcf, and in X-rays, are remarkably well correlated regarding their morphology and spatial extent. 

5_j Aims. In this work we present new results from a statistical analysis in order to specify and quantify these phenomenological relations. 

jy>j Methods. This is accomplished by investigating soft X-ray (0.3 -2.0keV) luminosities, FIR, radio continuum (1.4 GHz), Hor, 6-band, and UV 

C^ (1550A- I650A) luminosities for a sample of 23 edge-on late-type spiral galaxies. Typical star formation indicators, such as star formation 

k/ rates (SFRs), are determined and a statistical multi-parameter/frequency correlation analysis is carried out. 
• ,-H Results. We find strong linear correlations, covering at least two orders of magnitude, between star formation indicators and integrated 
ji^ (disk+halo) luminosities in all covered wavebands. In addition to the well established LFiR/ii.4GHz relation, we show new and highly signifi- 

?H ' cant linear dependencies between integrated soft X-ray luminosities and FIR, radio continuum. Ha, i?-band, and UV luminosities. Moreover, 
. . . , integrated soft X-ray luminosities correlate well with SFRs and the energy input into the ISM by SNe. The same holds if these quantities are 
plotted against soft halo X-ray luminosities. Only a weak correlation exists between the dust mass of a galaxy and the corresponding X-ray 
luminosity. Among soft X-ray luminosities, baryonic, and H i gas masses, no significant correlations are found. There seems to exist a critical 
input energy by SNe into the ISM or a SFR threshold for multi-phase halos to show up. It is still not clear whether this threshold is a physical 
one or represents an instrument dependent sensitivity limit. 

Conclusions. These findings strongly support our previous results that multi-phase gaseous galaxy halos in late-type spiral galaxies are created 
and maintained by outflowing gas produced in star formation processes in the disk plane. They conflict with the concept of halos being mainly 
due to infalling gas from the intergalactic medium. 

Key words. Galaxies: formation - Galaxies: halos - Galaxies: ISM - Galaxies: spiral - Galaxies: starburst - X-rays: galaxies 

1. Introduction sociated with diffuse ionized gas (DIG) in the halo. It is gen- 

j 1 1 1 erally accepted that the non-thermal radio continuum is a reli- 

In Paper! J TiillmannetaU |2006|) we established the morpho- ^^j^ indicator for cosmic ray (CR) halos, that extraplanar soft 

logical correlation and spatial coincidence between gaseous ha- x-rays are a good tracer of the hot ionized medium (HIM, 

los in late-type spiral galaxies, traced by their radio contin- ^ ^^ produced by supernova remnants (SNRs), in superbub- 

uum, Ha, and soft diffuse X-ray emission. We further showed ^j^^ ^^ -^ superwinds), and that the presence of extraplanar 

that the continuum emission radiated at UV wavelengths (~ pip is indicative of the so called disk-halo interaction (see 

210 nm) originated in the disks of the galaxies and is well as- \j^^„^^^ ^^ for ^ recent review). Therefore, we concluded 

. „ ; ; : , , . , ,^, ,>, ^t ^^^ • that multi-phase gaseous halos are created by star formation 

Based on observations obtamed with XMM-Newton, an ESA SCI- ^^„^ , , .,,.,, ^ , , 

■,, ■ , . , . u .■ J- .1 r J J u (SF) related processes in the disk plane. Our results presented 

ence mission with instruments and contributions directly funded by *^ *^ *^ 

ESA member states and NASA. i" P^P^"" ^ ^id not yield convincing evidence that galaxy halos 
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form b y infalling gas from an external reserv oi r as considered 
e.g. bylBenson et al.l(l2000l) . lToft et alJ(l2002l) or lPedersen elaP 
ll2005l) . 

Apart from these phenomenological analogies there also 
should exist strong statistical correlations between total 
(disk-Hhalo) soft X-ray, Ha (DIG), and radio continuum lumi- 
nosities. If this hypothesis is correct, these luminosities should 
also strongly correlate with other star formation i ndicators, 
such as the FIR, B-band, and UV luminosity (e.g., ICondonl 
Ll992; Read & Ponman, 2001), but also with direct tracers of 
massive st ar formation, such as star formation rates (SFRs, see 
[Kennicutr, 1998), supernova (SN) energy input rates normal- 
izedto the area of active st ar form ation £'°' (e.g.,Dahlem et al., 
Il995t iRossa & Dettmaii l2003[ simply called energy input 
rates, hereafter), and Lpi^/D-^^ ratios ( iLehnert & HeckmanI 
Il996l) . 

In the following we investigate the interrelations among 
star formation related key parameters and luminosities in or- 
der to establish significant (statistical) correlations of the up to 
now purely morphological relations. 

Integrated luminosities are used to test whether the diffuse 
emission from the disk and halo is coupled to SF, whereas the 
pure extraplanar soft X-ray emission is put into relation to the 
above mentioned parameters in order to show that multi-phase 
halos are indeed tightly correlated with SF processes in the 
disk. 

Finally, we discuss the question of the existence of an en- 
ergy threshold provided by stellar feedback processes which 
needs to be exceeded in order to create multi-phase halos. From 
the above mentioned correlations first empirical estimates of 
the suspected energy threshold are derived. 



overcome by selecting edge-on (/ > 70°) late-type spiral galax- 
ies from the literature for which the above mentioned luminosi- 
ties have been published and the SF related parameters could 
be calculated. In order to reach the largest possible coverage of 
the parameter space no constraints were imposed on the inte- 
grated luminosities, the S6Q/S 100 and Lpir/D^^ ratios, and the 
energy input rates E'°\ 

Galaxies are considered to have multi-phase gaseous halos, 
if they show kpc-sized extraplanar radio continuum, DIG, and 
X-ray emission which surrounds a large fraction of the disk 
plane. Extended X-sh aped DIG structures, as visible on Ho; 
images of NGC 1482 ("Strickl and et"ail |2QQ4a) or NGC 5775 
(Tiillmann et al., 2006 ), trace most likely a limb brightened 
outflow cone and are also classified to possess a halo. Single 
plumes or filamentary structures are not considered to consti- 
tute a halo. 

It should be kept in mind that the sample is limited in sen- 
sitivity and is not complete. 

T he enlarge d sampl e consists now of 23 galaxi es, ii ii ie frorn 
Tull mann etall ( l2006l) . seven from [Strickland eTaP (l2004al 
three of their targets are also in our sample and their results are 
consistent with ours), and another seven galaxies taken from 
the literature. The number of targets is sufliciently large to 
test the individual correlations by means of a statistical anal- 
ysis. All physical parameters adopted throughout this study are 
listed for the whole sample in Table 1 . 

Integrated X-ray luminosities for the first nine galaxies 
listed in Table 1 are from Paper I, the next seven are taken from 
Strickland et al. (2004a), whereas the last seven stem from 
other literature sources (see notes to table). 



2. Statistical analysis of key parameters 

2.1. Key parameters and sample selection 

It is straightforward to include additional quantities which trace 
massive stars in the disk and the energy they provide for the 
interstellar medium (ISM) and to search for trends among all 
those parameters. Such supplemental key quantities considered 
here are the B-band and UV luminosites as well as SFRs (or 
the FIR SFR per unit area, expressed by Lpir/D2^) and energy 
input rates by SNe, represented by E^'. 

We also include variables which are not directly coupled to 
SF, but trace different neutral phases of the ISM, such as the 
dust mass and the H i mass of a galaxy. 

Moreover, the extent of the SN-driven gas and therefore the 
size of the multi-phase halo also dep ends on the gravitational 
potential of a galaxy (e.g., iMac Low & FeiTara, 1999). As a 
first order approximation of this quantity the baryonic mass of a 
galaxy (stars and H i) determine d from the A'-band TuUy-Fisher 
relation (Bell & d e Jond. 120011) is also included. 

As was already pointed out in Paper I, our initial sample is 
biased against starburst galaxies, that is towards higher SFRs 
or energy input rates E^'. In order to investigate the possible 
coiTelations, non-starburst (but actively star forming) galaxies 
need to be included which have to cover the intermediate and 
lower energetic ends of these parameters. This obstacle was 



2.2. Statistical analysis 

We performed least-squares fitting (assuming Y = mX + b) 
and Spearman rank-order correlation analysis, to test the sig- 
nificance of the correlation between the investigated pairs of 
parameters. As these quantities, for example Lx and Lfir, usu- 
ally vary only within a certain physical range, e.g., set by the 
number of stars in the disk, these variables are most likely not 
normally distributed and their interrelation needs to be investi- 
gated by means of rank-order coiTelation analysis. 

The slopes m, their intercepts b and b' (for m - 1 .0), the 
reduced x^i the Spearman rank-order coiTelation coefficients 
r,, ?s, the levels of significance of r^, and the corresponding 
p-values of r^ are listed in Table |2l In the case that both quan- 
tities are not linearly correlated (null hypothesis), t^ follows a 
Student's f-distribution with n - 2 degrees of freedom. The null 
hypothesis turns out to be wrong if Ifsl > fn-2,i-Q', where a is 
the level of significance, which we adopt to be 0.01 . The quan- 
tile f2i,o.99 of th e f-distribution, as tab ulated in relevant statistics 
textbooks (^e.g.. lHartung et alll2002l) . is 2.518. For all fj-values 
which are larger, the null hypothesis needs to be discarded 
(with a 99% confidence level) which means that all those pairs 
of parameters are strongly coiTelated. 

Moreover, if the p-values listed in Table|2lfor the integrated 
luminosities are less than the adopted a-level, the null hypoth- 
esis needs to be discarded, too. In the present case only H i and 



R. Tiillmann et al.: The multi-phase gaseous halos of star forming late-type galaxies 3 

Table 2. Fit parameters and correlation coefficients. Listed are the Spearman rank-order correlation coefficient r,, its correspond- 
ing significance 4, and the p-value of r^, together with results derived from linear regression (Y - mX + b) for the pairs of 
parameters plotted in Figs.^|5] Uncertainties of the gradient m and the intercept b are given on a Icr level, b' is the intercept if a 
slope of m = 1 is assumed. 
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baryonic masses do not fulfill these criteria and are therefore 
assumed not to correlate with soft X-ray luminosities. 

The same analysis was carried out for the extraplanar soft 
X-ray emission (Lx.h) and the derived parameters are also listed 
in Table|5] All correlations and their corresponding least-square 
fits are shown in Fig.|3l 

3. Results and Discussion 

3.1. Correlations with integrated luminosities 

Based upon statistical analysis of the sample strong linear cor- 
relations (rs > 0.67) between star formation indicators and inte- 
grated (disk-Hhalo) multi-frequency luminosities could be es- 
tablished. In addition to the generally accepted LFiR/ii.4GHz 
relation (e.g., deJongetal., 1985; Condon, 1992, see also 
Table 2 of this work), highly significant linear dependencies 
between integrated soft X-ray luminosities (0.3-2.0keV) and 
integrated radio continuum (1.4 GHz), B-band, UV, Ha, and 
FIR luminosities have been found (see Fig.Q' Moreover, the 
integrated soft X-ray luminosities correlate well with star for- 
mation rates and SNe energy input rates (Figs.^andl^Jl. These 
correlations unambiguously confirm that the diffuse radiation 
of the X-ray emission and the DIG, as well as the CRs detected 
in the disks and halos are indeed coupled to SF 

Although Ha, B, and UV-band luminosities are generally a 
good tracer of the SF activity, these quantities can be seriously 
affected by absorption internal to the galaxy. The typical scatter 
of about 1 dex visible in our correlations might be indicative of 



this effect. However, absorption appears to be unlikely to lower 
the luminosities in such a way that the correlations are seriously 
weakened or even destroyed. 

The strongly discrepant data point visible in the Lx/^uv- 
diagram represents measurements for one of the most extreme 
starburst galaxies, M 82. In order for M 82 to follow our corre- 
lation, a much higher (lower) UV (X-ray) luminosity of about 
1 dex is required. However, independent observations car- 
ried out in the UV ("Cod e & Welchril98 2': 'Rifatto et al.', 'l995 ; 
Hoopes et al.. 2005) and at X -ray energies ( Fabbiano et al 



' The latter two relations are not shown here, as the Ho- and FIR 
luminosities are directly proportional to the corresponding SFRs (see 



lKenmcuttlfT99l . 



1992; Strickl and et all l2004al) yield very consistent results 
which imply that this "deviation" is an intrinsic feature of this 
outstanding galaxy. 

More than 60% of the UV emissio n of M82 measure d 
at about 1530A is extraplanar in origin ( iHoopes et al.l l2005h . 
These authors claim that extraplanar star formation or a com- 
bination of photoionization and shock ionization either appear 
unlikely or cannot account for the halo emission in the UV. 
Instead stellar continuum radiation from the disk plane scat- 
tered at dust particles in the halo shall be the main mechanism 
causing the UV emission. Given that M 82 is a low mass spi- 
ral galaxy (Mbar = 1.86 x IO'^Mq, Table 1) dust is assumed 
to be driven relatively easily by the superwind to high extra- 
planar distances. This independently supports the "scattering"- 
hypothesis and could also explain the lacking UV emission, 
provided the (Lx, Luv)-i"elation is vaUd. 

Alternatively, the high X-ray and low UV luminosity of 
M 82 can be explained by gas within the starburst region which 
is thermalized beyond the UV excitation limit. In such a ex- 
treme case in which the diffuse soft X-ray emission is domi- 
nating the energy output, we expect a deviation from the linear 
correlation curve during phases of strong starburst activity. 
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Fig. 1. Diagnostic diagrams confirming the strong correlations between integrated radio continuum, B-band, UV, and soft X-ray 
(0.3 -2.0keV) luminosities. X-ray luminosities also correlate well with the energy inp ut rate by SNe {Eh) and with FIR and Ha 
SFRs. Circles address our sample, squares represent data from Strickland et al. (2004b), and triangles denote data, collected from 
the literature (cf. Table. 1). Filled (open) symbols refer to galaxies with (un)detected multi-phase gaseous halos. Solid lines are 
best fits from linear regression while dashed lines indicate the trend expected for a relationship of unit slope, r^ is the Spearman 
rank-order correlation coefficient. 



Apparently only a weak correlation between the dust mass 
of a galaxy and the corresponding soft X-ray luminosity exists. 
Dust, as a byproduct of SF, is correlated with the FIR luminos- 
ity and should therefore also correlate (at least weakly) with 
other SF indicators, such as X-ray luminosities. 

The non-correlation between Lx and Mh i is also not sur- 
prising as the neutral gas component of the ISM is not directly 
involved in SF. 

Surprisingly only at a first glance, a correlation between 
galaxy potentia l and X-ray lum inosities is not found. As the 
baryonic mass (iBell & de Jona 12001.) is calculated from ro- 
tational velocities of the stars and the H i-gas in a galaxy, this 



parameter traces the dynamics of the system and is therefore in- 
sensitive to SF. Nevertheless, this finding bears two important 
implications. Firstly, the dynamics of the gas is apparently not 
related to the diffuse emission produced in SF events. Secondly, 
for a more accurate estimate of the baryonic mass a better mass- 
related approximation needs to be found. This relation should 
also include dark matter halos in order to study the effects this 
additional component might have on the established relations. 
If there is still no correlation with baryonic mass, then the tar- 
gets of our sample do not foll ow the Lx ~ v^„, relation derived 
from the infall simulations bv lToft et alJ 1120021) . 
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Fig. 2. Functional dependence between soft X-ray luminosities and Lfir/ZJjs (upper left), the dust mass of a galaxy (upper right), 
its H I mass (lower left), and the baryonic mass (lower right). 
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Fig. 3. Functional dependence between Lx,h, the soft (0.3 - 2.0 keV) X-ray luminosity of the halo, and Lyi^jD\^ (upper left), ^a 
(upper right), S FR^^^r (lower left), and S F/?h» (lower right). Solid symbols represent our sample, open symbols denote data taken 
from aStrickland et al. ( 2004a) . Solid and dashed hnes have the same meaning as described in Fig.Q] 
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It should be pointed out that Mbai represents the gravita- 
tional potential of a galaxy and should therefore be sensitive 
to the overall extent of the multi-phase halo. However, it is not 
sensitive to the existence/non-existence of such halos. In this 
regard the SFR and the energy input rate by SNe into the ISM 
turn out to be the most important parameters. 

3.2. Correlations with extraplanar luminosities 

In order to establish that multi-phase halos are indeed a product 
of SF, the pure extraplanar emission in the radio continuum, 
in Ha, and at soft X-ray energies also needs to coiTelate with 
SFRs and the energy input rate by SNe. 

At present this analysis can best be done in the X-ray 
regime. Other wavelengths, such as Ha and the radio contin- 
uum, are not usable, because the integral disk and halo lumi- 
nosities do not allow to reliably estimate the fraction of the lu- 
minosity originating in th e halo. For the DIG, this fraction can 
vary from 12% to 60% feossa & Dettmaii I2OO0I) and would 
therefore lead to unacceptable uncertainties. 

From the original sample, 13 galaxies are suitabl e for 
the "halo"-sample (those from IStrickland et al.l ll2004al) and 
[riillmann et al. (2006)). The remaining galaxies either do not 
show extraplanar X-ray emission or the published data did not 
allow a disentangling of disk and halo luminosities. 

It should be pointed out that inclination issues, regarding 
the reliability of detecting multi-phase halos, are negligible. 
Either because extraplanar gas was detected in case of incli- 
nation angles < 80° (see NGC 1482, Strickland et al., 2004a), 
or because the inclination-corrected extent of the "extraplanar" 
gas was indeed not large enough to be seen above the disk plane 
(see NGC 3877, Tullmann et al., 2006). Most importantly, Lx 
and Lx,h are not linearly correlated with the inclination angle 
i. This is proved by low Spearman rank-order coiTelation coef- 
ficients of rs - -0.26 and r^ - 0.23, respectively, again on a 
99% confidence level. 

In Fig. |3]we plot the soft (0.3-2.0keV) halo luminosity 
(Lx,h) as a function of the different SF indicators (which are 
SFRpij^, SFRiia, Lpyr/D^^, and Ea)- All pairs of parameters 
show strong correlations (r^ > 0.67) on a 99% confidence 
level, except Lx,h vs. Lps^/D^^ (rj - 0.61). A possible expla- 
nation could be that Lpiji/D^^, a measure of the SF rate inten- 
sity per unit disk area, underestimates the real SF intensity. In 
disk galaxies SF usually occurs within the Hor-disk and is not 
that widespread as the distribution of the old stellar population 
(measured by D25) implies. 

It turns out that the FIR and the Ha SFRs are the most im- 
portant parameters which decide on the formation of X-ray ha- 
los. Among these two quantities, 5F7?fir seems to be the more 
reliable parameter, because it is very little aff'ected by dust ab- 
sorption and is dominated by disk emission. On the other hand, 
the Ha emisson, and therefore S FRua, is a more sensitive mea- 
sure of how much gas is converted into massive stars. However, 
this quantity is significantly affected by extinction and in some 
cases most likely dominated by extraplanar emission. 

Given the strong correlations between integrated X-ray, 
Ha, and radio continuum luminosities as well as between halo 



X-ray luminosities and SFRs or Ep,, it is plausible, to adopt 
similar relations for the extraplanar DIG and the CR halos. 
Therefore, we expect the SFR and the energy input rate by SNe 
to be the main driving-agents for multi-phase galaxy halos. Our 
statistical results are also fully consistent with the morpholog- 
ical and spatial coincidences presented in Paper I. For consis- 
tency reasons it is certainly worthwhile to re-examine existing 
radio and DIG data to allow for a clear separation of disk and 
halo luminosities. 



3.3. Constraining the star formation threshold 

In order to achieve a more comprehensive picture on the evolu- 
tion of galaxy halos in different wave bands, we need to answer 
the question of the existence of a threshold SFR (or equiva- 
lently the SN energy input rate Ea) required to create multi- 
phase gaseous halos. The existence of such a threshold can be 
motivated as follows: It is well known that hot ionized gas is 
driven off-plane via superbubbles produced by SNe. From a 
theoretical study of breakout conditions of such superbubbles 
f Mac Low & Mc Cravll 19881) we know that the mechanical lu- 
minosity of SNe is a key quantity which determines whether 
a superbubble will breakout or collapse. In other words, the 
overpressured X-ray emitting gas will leave the disk and form 
a halo if the energy provided by SNe is above a certain limit. 

On the other hand, it can be argued that channels, through 
which the hot gas reaches the halo, were already blown into the 
ISM during previous periods of intense star formation. Hence, 
a substantially lower SFR is needed to elevate the gas into the 
halo, making it even more difficult to establish a threshold by 
means of observations. 

From Fig.[nthe existence of a critical SFR (or energy input 
rate) is hard to constrain, as the lower energy end of our coiTe- 
lations is statistically not well covered. It appears, however, that 
galaxies without halos (open symbols) also follow a linear rela- 
tion and that the region dividing galaxies with halos from those 
without is relatively narrow. For multi-phase halos to evolve, 
the data imply that a critical threshold of S F/?fir > 1 Mo/yr, 
SFRua > 0.1 Mo/yr or ^a > 1 x 10"^ erg s"' cm"^ needs to 
be exceeded. Interestingly, DIG seems to coexist with other gas 
components fo r Ha energy fluxes large r than (3.2 + 0.5) x 10"^" 
erg s ' kpc"2 (JRossa & Dettmaill2003l) . 

Unfortunately, with the present data a detection threshold, 
introducing systematic uncertainties by the limited sensitivity 
of the instrument, cannot be a priori ruled out. In order to dis- 
tinguish between a significant physical threshold and a sensi- 
tivity limit as well as to reach a better understanding on the 
formation of gaseous galaxy halos, deep observations of galax- 
ies with no or only little halo emission are required. Therefore, 
future studies should aim at increasing the sample to test the 
low energy end of our coiTelations and to investigate the above 
relations for galaxies of different Hubble-types, such as Sa and 
Irr, as the different SF histories of these galaxies might lead to 
different relations and thresholds. 
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4. Summary and Conclusions 

With a sample of 23 actively star forming galaxies studied in 
the X-ray regime and by implementing additional wave bands 
(Ha and UV), we found remarkably strong linear coiTelations 
between integrated 1 .4GHz radio continuum, FIR, Ha, B-band, 
UV, and soft X-ray luminosities. Strong correlations also exist 
if soft X-ray luminosities are plotted against SFRs, Lfir/£'25' 
or the energy input rate by SNe per unit area, expressed by Ep^. 
Integrated X-ray luminosities neither correlate with the H i nor 
with the baryonic mass of a galaxy. 

Strong correlations are also found if the diffuse soft X-ray 
luminosity of the halo is plotted against the FIR and Ha SFR 
or the SN energy input rate Ej^. These quantities are considerd 
to be the most important parameters for the creation of multi- 
phase halos. 

If a critical energy threshold exists and an instrumental 
detection bias is negligible, the present data suggest that a 
threshold of 5 F/?FiR > l-QM^Iyr, SFR^a > O.lMo/yror 
Ea > 1 X 10""' erg s"' cm"^ needs to be exceeded in order 
to create multi-phase galaxy halos. 

Our results are in agreement with previous findings of a 
morphological and spatial coincidence of gaseous multi-phase 
halos (see Paper I). They clearly imply that multi -phase halos 
are the con sequence of stellar feedback processes i n the disk 
plane (e.g.. Ide Avillez & Breitsch werdt*. '2004*. '200^) but con- 
flict with the concept of halos being due to infallin g gas from 
the intergalactic medium (e .g.. Benson et all l200(i iToft et all 
l2nn2HPedersen et allEool . 
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